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ABSTRACT 

We search the BOOMERanG maps of the anisotropy of the Cosmic Microwave Background (CMB) for 
deviations from gaussianity. In this paper we focus on analysis techniques in pixel-space, and compute 
skewness, kurtosis and Minkowski functionals for the BOOMERanG maps and for gaussian simulations 
of the CMB sky. We do not find any significant deviation from gaussianity in the high galactic latitude 
section of the 150 GHz map. We do find deviations from gaussianity at lower latitudes and at 410 GHz, 
and we ascribe them to Galactic dust contamination. Using non-gaussian simulations of instrumental 
systematic effects, of foregrounds, and of sample non-gaussian cosmological models, we set upper limits 
to the non-gaussian component of the temperature field in the BOOMERanG maps. For fluctuations 
distributed as a 1 DOF x 2 mixed to the main gaussian component our upper limits are in the few % 
range. 

Subject headings: Cosmic Microwave Background 



1. INTRODUCTION 

In most inflationary scenarios the primordial density 
field is expected to be gaussian (see e.g. Peebles 1999). In 
contrast, structure formation scenarios based on topolog- 
ical defects (e.g. Avelino et al. 1998) or less general in- 
flationary models (Linde et al. 1997; Martin et al. 2000; 
Contaldi et al. 1999) predict a non-gaussian density field. 
Thus, measurement of the statistical nature of the 
CMB anisotropies can distinguish between these scenar- 
ios (see e.g. Fischler et al. 1985). 

The presence of noise in the measurements combined 
with the limited coverage of the present observations 
can mask cosmological non-gaussian features. Moreover, 
the presence of systematic effects can produce subtle in- 
strumental non-gaussian features in intrinsically gaussian 
anisotropy maps. 

Efforts to identify non-gaussianities in the CMB 
have been extensively carried out for the COBE 
data (see e.g. Ferreira et al. 1998; Banday et al. 1999; 
Komatsu et al. 2001; Bromley and Tegmark 2000) and no 
significant detection of cosmological non-gaussianity has 
been reported (but see Magueijo 2000). However, the sen- 



sitivity of COBE-DMR to the expected levels of cosmo- 
logical non-gaussianity from rare highly non-linear events 
like topological defects (see e.g. Durrer and Zhou 1996) 
is not very high, due to the large field of view which 
smears out the effects of small scale features. While the 
present observations on the power spectrum are incon- 
sistent with a structure formation scenario solely based 
on defects (see e.g. Durrer et al. 1998), a mixed infla- 
tion+defects model is still compatible with the data 
(Bouchet et al. 2001). 

Furthermore, the low signal-to-noise ratio and coarse 
resolution in the COBE maps also makes it difficult 
to detect primordial non-gaussianity; the central-limit 
theorem states that the sum of various instrumental 
effects will make the distribution tend to a gaussian 
(Novikov et al. 2000). 

Analyses to date of higher angular resolution exper- 
iments, like QMASK (Park et al. 2001) and MAXIMA 
(Santos et al. 2001; Wu et al. 2001), show full consistency 
with gaussianity. 

In this paper we focus on the maps produced by 
the BOOMERanG experiment (de Bernardis et al 2000; 
Netterfield et al. 2001). Due to its wide sky and frequency 



1 



2 



coverage, BOOMERanG is ideally suited to carry out an 
accurate analysis of the possible systematic effects present 
in the detected signal. We analyze them using a Monte- 
Carlo approach, in order to set quantitative upper lim- 
its for the level of primordial non-gaussian fluctuations 
present in the maps. 

The data we use are presented in section 2. We use 
five pixel-space estimators of departures from gaussian- 
ity: the skewness and kurtosis of the CMB temperature 
distribution, and the three Minkowski functionals: area, 
length and genus (Minkowski 1903). In section 3 we apply 
these estimators to the BOOMERanG maps and to gaus- 
sian Monte-Carlo simulations of the CMB sky observed by 
BOOMERanG. In section 4 we analyze template maps 
of other non-gaussian signals which could in principle be 
present in the BOOMERanG maps; we study how sensi- 
tive the pixel-space techniques are in detecting these ef- 
fects, and we estimate how much these affect the measure- 
ments of cosmological non-gaussian signals. In section 5 
we compare our measurements to the cosmological expec- 
tations for a few sample scenarios. 

The methods presented here are especially useful in 
the detection of highly non-linear features (as expected 
in topological defects theories or for Galactic foregrounds) 
(Phillips and Kogut 2001). One may alternatively use a 
wavelets approach (Hivon et al. 2002), or a bispectrum 
approach (Contaldi et al. 2002). 

2. THE DATA 

We use the BOOMERanG data obtained in the Long 
Duration Balloon (LDB) flight in 1998. For this paper 
we try to be as conservative as possible, so we concen- 
trate on the center of the map, 1.19% of the full sky 
(70° <RA< 105°, -55° < S < -35°), where the sky cov- 
erage is approximately uniform and the integration time 
per pixel is maximum. The galactic latitude of this region 
spans from —42° to —13°. We also consider for compar- 
ison a lower latitude region with similar size but extend- 
ing to b < —5°. We concentrate on the best 150 GHz 
channel, B150A, for CMB signals, and on the best 410 
GHz channel, B410B1, for monitoring the dust foreground. 
We use only data collected in ldps scans. Since pixel- 
space estimators are sensitive to the size of the pixel, we 
have used 28', 14', and 7' pixels (HEALPix pixelisation 
(Gorski et al. 1998)). In fact, as the pixel size is made 
smaller, smaller size structures become visible, but the 
signal to noise ratio (SNR) per pixel decreases. In the 
reference region, the noise per pixel at 150 GHz is of the 
order of 30, 50, 90 [iK for 28', 14', and 7'pixels respec- 
tively. For this reason it is important to investigate the 
best tradeoff between angular resolution and SNR. The 
best fit dipole anisotropy component is removed from the 
time ordered data (TOD) before further processing. TOD 
are deconvolved from the frequency response of the in- 
strument and high-pass filtered with a sharp filter, with a 
cut-on at 70 mHz. We use the improved pointing solution 
of (Netterfield et al. 2001), featuring an equivalent beam 
FWHM of 12.7' for the B150A channel. The time ordered 
data obtained in this way are converted into a map using 
the fast iterative generalized least squares map making of 
(Natoli et al. 2001). The time-domain high-pass filtering 
procedure removes most of the 1/f noise and most of the 



scan-synchronous noise (SSN, see section 4) , but also intro- 
duces an important anisotropy in the large-scale content 
of the map: large scale components of the sky tempera- 
ture distribution are filtered out in directions parallel to 
the scan, while being maintained in directions orthogonal 
to the scan. This is not completely recovered by the map 
making. Since the scan direction is not constant during 
the flight, it is difficult to model analytically the effect of 
such a filter. For this reason all the subsequent analysis is 
performed numerically, using Monte-Carlo techniques. 

3. ANALYSIS OF THE MAPS 

The normalized skewness S3 and kurtosis S4 are com- 
puted from the 2 nd , 3 rd and 4 th moments a 2 = Y,ii T i ~ 
(T)f/(N 1), M3 = - {T)f/N and ^ = £ 4 (T, - 

{T)) 4 /N of the distribution of the T t : S 3 = H 3 /a 3 , 
S 4 = M4/ff 4 -3. We find S 3 . m = -0.07 and 5 4 , m = -0.10. 
These values are plotted as vertical lines in fig.l. As 
a comparison, the 410 GHz channel in the same region 
produces S 3 . m = 0.32 and S^ m = 1.01. Due to the 
presence of sample variance, as well as noise, correlations 
and anisotropic filtering, these values are not expected to 
be exactly zero, even for a perfectly gaussian sky. For 
this reason we compare the measured data to Monte- 
Carlo simulations. These are based on gaussian maps of 
the CMB sky extracted from a parent distribution with 
the best fit power spectrum measured by BOOMERanG 
(Netterfield et al. 2001). These maps are sampled along 
the BOOMERanG scans, simulated noise with the statis- 
tical properties of the considered channel is added, and 
the resulting TOD is filtered in the same way as the 
BOOMERanG data. The SSN is not included in the sim- 
ulation, but its contribution is negligible (see section 4). 
We then obtain simulated observed maps using the same 
algorithm we have used for the actual BOOMERanG sig- 
nals. From each simulated map we compute 53 and S4, 
and plot their distribution in fig.l. 




FIG. 1. Measured values of the normalized skewness (top) and 
kurtosis (bottom) in the nominal box (left) and in the low lat- 
itudes box (right) at 150 GHz are shown as vertical lines. The 
histograms are computed from detailed Monte-Carlo simula- 
tions of the measurements. The simulations were obtained from 
a combination of a gaussian sky with the best fit power spec- 
trum measured in B98 and a simulation of detector noise. The 
different width of the distributions is due to uneven coverage 
and noise per pixel. 

The visual comparison of the actual measurement to the 
simulations suggests consistency with a gaussian sky at 150 
GHz. From the distributions we compute the single-sided 
probabilities P(S 3 > S 3>m ) = 79% and P(S 4 > S A . m ) = 
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73% for the 150 GHz channel, while for the 410 GHz chan- 
nel P(S 3 > S 3 , m ) = 2% and P(S 4 > S 4 , m ) = 1% • In the 
region at lower Galactic latitudes P(S 3 > S 3m ) = 73% 
and P(S 4 > S i>m ) = 2% at 150 GHz, and P(S 3 > S 3 , m ) < 
10~ 8 and P(S 4 > S A , m ) < 10~ s at 410 GHz. All the re- 
sults above are for a 28' pixel size. Similar results are 
obtained with 14' and 7' pixels. 

For a given threshold v we compute the surface densities 
of three Minkowski functionals Vj following the contours 
dQ of the excursion sets Q = Q{v) = {Ti : {Ti - (T))/a > 
v} (see e.g. Gott et al. 1990): 

V ° = ^J Q dA > V ^4l/ a / S ' V ^^AS dQ MS W 

where dA and ds are the differential elements of Q and 
dQ respectively, and k is the geodesic curvature of ds. 
The values of the v^ m measured from the 150 GHz map 
are plotted as diamonds in fig. 2. We show the results for 
28' pixels, since the CMB signal is larger than the noise in 
this case. Again, even for a gaussian sky, the presence of 
noise, correlations and anisotropic filtering does not allow 
us to compare the measurements to the analytic expres- 
sions of Vj which can be computed for gaussian isotropic 
distributed data. So we compute the distributions of the 
vj from the simulated maps. We plot the 95% confidence 
band of the Vj vs threshold v in fig. 2. 
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FIG. 2. Measured values of the area (top), length (middle) and 
genus (bottom) Minkowski functionals in the nominal box at 
150 GHz are plotted as diamonds vs the normalized thresh- 
old. The solid lines represent the 95% confidence band of 
the distributions of the same quantities as computed from de- 
tailed Monte-Carlo simulations of the measurements, obtained 
assuming a purely gaussian sky with the best fit power spec- 
trum measured in B98, and 5a realistic simulation of detector 
noise. Dashed lines are 5analytical predictions for the gaussian 
isotropic case. A 28' pixelization has been used. 

The consistency of the measured Vj iTn with the gaussian 
sky hypothesis is evident. At this pixelization, the Vj. m 



are very similar to the ones computed from CMB only 
simulations. At smaller pixelizations the effect of noise is 
to increase Vi and v 3 with respect to the CMB-only case, 
reducing the cosmological significance of the results. 

To compare numerically the measured Minkowski func- 
tionals to the simulations, we define 

X? = E [<» - <<"»>] • [< m M,< m (M)] • 

■[<>)-<<"(/*)>] 

where c [w| lm (^), w| mi (^)] is the covariance matrix. From 
the measured data (with 28' pixels) we have Xo = 9, 
X\ = 7, xl = 15, and P( X 2 > xl) = 76%, P( X 2 > 
Xi) = 91%, P(x 2 > xl) = 31%: again an excellent agree- 
ment with the gaussian sky hypothesis. As a compari- 
son, for the 410 GHz data we have Xa — 23, xi = 8, 
X 2 2 = 16, P( X 2 > Xl) = 0.1%, P( X 2 > Xl) = 84%, 
P(x 2 > Xi) — 26%. In the low Galactic latitudes box 
the Gaussian hypothesis is rejected for all the functionals, 
for both the 150 GHz channel and the 410 GHz channel. 
Similar numerical results are obtained for 14' and 7' pixels; 
we note, however, that for the finest resolution the results 
are noise dominated rather than CMB dominated. 

4. ESTIMATES OF SYSTEMATIC EFFECTS 

Though we do not detect any evidence for non- 
gaussianity in the maps, we set upper limits on possible 
sources of non-cosmological non-gaussianity in order to as- 
sess the potential for further analysis using more data. The 
raw BOOMERanG TOD are affected by a large scale SSN, 
typically of the order of 20 in the 150 GHz channels 

and of 400 fiK/° (CMB temperature units) in the 410 
GHz channels. The delay between the SSN and the az- 
imuth observed suggests an instrumental, thermal origin 
of the SSN. We have subtracted the best fit CMB dipole 
from the 150 GHz TOD and low passed the TOD to keep 
only the large scale SSN. We have then built a map from 
this signal with the usual map making algorithm, includ- 
ing the high pass filter, which strongly reduces the SSN 
amplitude (see fig. 3). Despite the relatively large signal in 
the TOD, the anisotropics induced in the filtered map are 
of the order of few /iK in the central box considered for 
this analysis. 

These are very small with respect to the detected 
anisotropies, but evidently non-gaussian. We have used 
this map as a template of a systematic effect. We built 
simulations of contaminated data by adding this map to 
the gaussian simulations of the CMB used in section 3. 
Then we have applied our gaussianity estimators to the 
simulated maps. The result is that the distributions of 
the estimators are not affected by the presence of SSN at 
the level we had in the 150 GHz channel analyzed here. 
We compare the x 2 of the gaussian-only simulations de- 
fined above to the x 2 of the contaminated simulations, for 
different levels of the SSN map. We find that the mean 
square amplitude of the SSN fluctuations should be of the 
order of 35 to 70% of the mean square CMB fluctuations 
(depending on the estimator) before becoming detectable 
at 95% C.L.. This means that the SSN contribution would 
have to be at least 7 times larger than that actually ob- 
served to affect the results reported in section 3. 
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FIG. 3. Estimated systematic effects at 150 GHz in the ob- 
served sky region. The top panel shows the effect of Scan Syn- 
chronous Noise (SSN) when projected on the sky using the same 
map-making procedure used for the actual data. The middle 
panel represents the extrapolation of interstellar dust emission 
from the FSD maps. The bottom panel shows the expected sig- 
nal from extragalactic point-sources extrapolated at 150 GHz 
using the WOMBAT catalogue. 28' pixels have been used. 

The second source of a systematic effect we have consid- 
ered is thermal emission from interstellar dust. We have 
used the FSD maps (Finkbeiner et al. 1999) filtered as in 
(Masi et al. 2001) to produce a template 150 GHz dust 
map. This is strongly non-gaussian, as visible in fig. 3, but 
its level is very small with respect to the general level of 
CMB anisotropy. We have added this map to the gaus- 
sian CMB simulations, and applied our gaussianity tests. 
Again, the \ 2 and the distributions of the contaminated 
simulations are very close to those of the gaussian-only 
simulations. Interstellar dust emission would have to be 3 
times larger than the FSD prediction to contaminate the 



results reported in section 3 at 95% C.L.. 

The third source of a systematic effect we have consid- 
ered is emission from extragalactic sources present in the 
observed region. We have used the WOMBAT catalogue 
(2001) to extrapolate to 150 GHz the flux of known radio 
sources in the surveyed area. From the catalogue obtained 
in this way we create a map convolving the point-source 
flux with the angular response of our telescope. The map 
(fig. 3) — which is evidently non-gaussian — is added to 
the gaussian CMB simulations and the analysis with our 
estimators is carried out. Yet again, the y 2 an d the distri- 
butions of the contaminated simulations are very close to 
those of the gaussian-only simulations. Fluxes from extra- 
galactic sources would have to be 6 times larger than the 
WOMBAT predictions to contaminate the results reported 
in section 3 at 95% C.L.. 

5. COSMOLOGICAL SOURCES OF NON GAUSSIANITY 

There are many models of non-gaussian CMB fluctu- 
ations. We have considered a few which can be easily 
compared to our results. As a first concrete example, we 
have used Linde's model (Linde et al. 1997). In this case 
CMB fluctuations are distributed as a y 2 with 1 DOF. This 
hypothesis is strongly rejected by all of our estimators. 
In other models there is a mixture of gaussian and non- 
gaussian CMB anisotropies (e.g. (Bouchet et al. 2001)). 
For example, Lyth and Wands ((2001)), consider a mix- 
ture of 1 DOF x 2 distributed fluctuations with Gaussian 
CMB anisotropies. From Monte Carlo simulations we find 
that the 95% confidence upper limits to the rms ampli- 
tude of the non-Gaussian signal are 5%, 6% and 6% of the 
Gaussian one (using Area, Length and Genus Minkowski 
estimators respectively, with 28' pixels); the upper limits 
become 2%, 5% and 8% respectively, with 7' pixels. The 
upper limit we get is not very sensitive to the pixelization 
because the signal to noise ratio increases by increasing 
the pixel size, but the angular resolution decreases, and 
there are less spots on which the functionals can be built. 
The two effects roughly compensate each other. 

We have also simulated CMB maps with a speci- 
fied (equal-£ reduced) bispectrum Bi = AC^ 2 (see 
(Contaldi and Magueijo 2001)). For C e we used the B98 
best fit power spectrum. For the non-Gaussianity estima- 
tors considered in this paper, we detected no deviation 
from Gaussianity significant at the 95% CL over the range 
of amplitudes A from to 10. 

6. CONCLUSIONS 

The pixel-based techniques discussed here and applied 
to the central region observed by BOOMERanG at 150 
GHz confirm the Gaussianity of the detected CMB fluctu- 
ations, and exclude with high confidence \ 2 distributed 
CMB temperature fluctuations. These techniques also 
detect non-Gaussian fluctuations due to interstellar dust 
at high frequencies and at lower latitudes. Subdominant 
non-Gaussian fluctuations mixed to the Gaussian ones are 
strongly excluded for some models, while are not detected 
at all in other cases. 
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